typeset using JPSJ.sty <ver.l.Ob> 



Charge and Orbital Ordering in the Triangular-Lattice £ 2£ rOrbital 

System in One Dimension: 



A Possible Ground State of BL V s O 



16 



o 
o 

(N 

> 

o 
O 
o 



c3 



o 
o 



> 

OS 



o 

03 



c 

o 
o 



X 
S3 



Yoshinori Shibata 1 * and Yukinori Ohta 1,2 ** 

1 Graduate School of Science and Technology, Chiba University, Chiba 263-8522 
2 Department of Physics, Chiba University, Chiba 263-8522 

(Received September 25, 2001) 

We consider the possible charge and orbital ordering in a Hollandite compound Bi^VsOie, 
which is a new one-dimensional triangular-lattice £2 9 -orbital system. Using the strong-coupling 
perturbation theory, we derive the effective spin-orbit Hamiltonian in the approximation ne- 
glecting the small off-diagonal hopping parameters or orbital fluctuation, whereby we obtain 
the spin Hamiltonians in the partial space of each orbital-ordering pattern. We then apply an 
exact-diagonalization technique on small clusters to these spin Hamiltonians and calculate the 
ground-state phase diagram. We find that a variety of orbital-ordering patterns appear in the 
parameter space, which include the state characterized by the partial singlet formation consistent 
with recent NMR experiment. 

KEYWORDS: charge ordering, orbital ordering, vanadate, Hollandite, Bi^VgOie, perturbation theory 



, §1. Introduction 



Orbital physics in transition-metal oxides, such as 
the orbital ordering (00) of the e g -spins of Mn-oxides, 
has attracted much attention in the research field of 
strongly correlated electron systems. 1 ' Recently, Kato 
et al. 2 ^ measured the magnetic susceptibility and resis- 
tivity of Bi^VsOie with 1.72 < x < 1.8 and found a 
metal-insulator transition (MIT) at T < 80 K. The crys- 
tal structure of this compound belongs to a group of 
Hollandite- type phases and has a VsOi6 framework com- 
posed of double strings of edge-shared V06 octahedra as 
shown in Fig. 1. Kato et al. 2 ^ suggested a mechanism 
of the transition that the charge ordering (CO) between 
V 3+ and V 4+ accompanied by an OO occurs with MIT. 

In the viewpoint of orbital physics, this system may 
be regarded as a one-dimensional (ID) version of L1VO2 
known as a possible 00 system of i2g-orbitals on the 
2D triangular lattice of S = 1 spins, 3,4 ' although our 
system Bi^VsOie has the average valence of \A 3+1 ' 3 ' + 
at x = 16/9 and thus is in the mixed valent state of 
V 3+ . V 4+ = 3d 2 . M i = 2 . 1 2) xhe centra i issue 

in the present system is therefore the mechanism of the 
MIT concerning how the highly frustrated spin, charge, 
and orbital degrees of freedom at high temperatures are 
relaxed by lowering temperatures and what type of the 
ground state is realized at zero temperature. 

To consider this issue, we study in this paper the possi- 
ble ground state of this system. We thereby assume that 
the V-ions, which have the triply degenerate ^g-orbitals 
with intra- and inter-orbital Coulomb interactions as well 
as Hund's rule coupling, form the ID triangular lattice 
and that the intersite Coulomb repulsions are strong 
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Fig. 1. Schematic representation of the double string of edge- 
shared VO6 octahedra in Bi^VsOig. 



enough for the electrons to be localized to form a spatial 
CO pattern. We then use the strong-coupling perturba- 
tion theory 5 ' to derive the effective spin-orbit Hamilto- 
nian, which turns out to be block-diagonal with vanish- 
ing orbital-off-diagonal elements within the approxima- 
tion used. We employ a numerical exact-diagonalization 
technique on small clusters of this Hamiltonian to con- 
sider the possible 00 spatial patterns in the ground 
state. We find that a variety of the ground-state phases 
appear in the parameter space. We also discuss the spin 
degrees of freedom of the obtained 00 patterns. We ar- 
gue that the state characterized by the singlet formation 
of two S = 1 spins on the neighboring V 3+ -ions with re- 
maining nearly-free 5 = 1/2 spins on the V 4+ -ions might 
be relevant in the present material. 

Although the experimental data on this new system 
Bi^VsOie are quite limited at present, we hope that our 
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first theoretical study on its charge, orbital, and spin 
degrees of freedom would stimulate further experimental 
studies of this intriguing material. 

§2. Model 

Our starting high-energy Hamiltonian is of the follow- 
ing form: 

H = H + H t 

(ij) iaa' ,OL^0 

i,a.^{3 ia. 

H t = - U a ,jp(4 aa c jf}a + H.c.) 

(ia ,jP).a 

where V is the intersite Coulomb repulsion, Jh is the 
Hund's rule coupling, and U and V are the intra- 
and inter-orbital Coulomb repulsions, respectively. (• • •) 
stands for the nearest-neighbor pair of sites. We neglect 
the crystal-held splittings among the ^g-orbitals for sim- 
plicity because the gain in kinetic energy may be much 
larger than the splittings as has been assumed in ref. 3 - 1 
We also assume the relation U' — U — 2Jh, which is 
valid in the atomic limit. 6 ) ti a .jp is the hopping param- 
eter between the orbital a on site i and orbital (3 on 
site j where a, (3 £ {d xy , d yz , d zx } in the coordinate sys- 
tem shown in Fig. 1. We retain the direct hoppings be- 
tween the i2g-orbitals on the V-ions because the indirect 
hoppings via the O-ions are rather small. ^ Independent 
nearest-neighbor ti a ^p parameters are t a , £fc, and t c as 
shown in Fig. 2. c| Q(T (ci QCT ) is the electron creation (an- 




Fig. 2. Schematic representation of the t2 9 -orbitals on the ID 
triangular lattice. Two of the four lobes for each of the three 
t2 9 -orbitals are drawn. 



nihilation) operator at site i, orbital a, and spin <r, and 
ni a a — c l a a c iaa is the number operator. We also define 
n ia = riiai + n lai with a =|, J. and = ^ a n la . 

§3. Charge ordering 

First, we consider the CO in the ground state. We 
assume the limit of no doubly occupied orbitals. Then, 
if we assume that the inter-site Coulomb repulsions V (as 
well as the inter-orbital Coulomb repulsion U') is much 
larger than the hopping parameters, we find the ground 
state of the system to be charge ordered. It is readily 
noticed that the lowest-energy CO state has the ordered 
pattern like • • • V 3+ V 3+ V 4+ • • • as shown in Fig. 3. At 
H t = 0, this state has the ground-state energy 

N 

E = —(4U'-4J H +32V) 
6 

when there are no doubly occupied orbitals, whereas the 
states containing V 2+ or V 5+ have the ground-state en- 
ergy 

N 

E = tt(6{7' -6 J H + 28V) 
o 

for ...V 2+ V 4 +V 4+ ---, and 

E = ^(8U' -8J H + 2AV) 
b 

for ■ ■ • V 2+ V 3+ V 5+ • ■ ■. Since the states containing V 2+ 
or V 5+ are highly unrealistic, we should have the condi- 
tion 

W-2V- J H > 
for the presence of CO of only V 3+ and V 4+ ions. 




Fig. 3. Schematic representation of the CO pattern in the ID 
triangular lattice. 



The ground state of the unperturbed (H t = 0) Hamil- 
tonian is then of the degeneracy 

M = 3 N ■ 3 2N/3 ■ 2 N ' 3 

because there are 3^ choices of orbitals, 3 2JV / 3 choices of 
5=1 spins, and 2 N ' 3 choices of S = 1/2 spins, where N 
is the number of sites in the system. The degeneracy is 
lifted by the perturbation of small hopping parameters 
tia.jp as shown in the next section. 

§4. Effective Hamiltonian 

We carry out the second-order perturbation calcula- 
tion with respect to Uajp assuming that ti a ,jS3 is much 
smaller than U, U', and V. The degeneracy is lifted by 
the perturbation and the effective spin-orbit Hamiltonian 
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is obtained as follows: 

ffeff = go -EiM)E <Ml TT lM,) ^ 

where \fi) and (/i,/x' = 1, ■■■ ) M) are the M inde- 
pendent eigenvectors of the ground state of Hq and \n) 
are the n-th excited states of Hq. Eq and E n are the 
corresponding eigenenegies of Hq. 

To find the expressions for the effective Hamiltonian, 
let us first prepare the degenerate eigenstates of Hq: 

iw=^c(M)nk^,»„«, d ) M 

where Si and sf are the spin quantum numbers at site-i, 
rii is the number of electrons at site-?, nf is the number 
of doubly occupied orbitals at site-i, and {ai\ is the set 
of occupied orbitals. {q} represents a set of all these qan- 
tum numbers. The basis states \si, sf , rii, nf){ a .} may be 
of the following form: 

(i) For the d 1 -site, we have 

11/2,1/2,1,0)* = | 1) a 
1 1/2, —1/2, 1, 0) Q = | |) Q 
with 5=1/2 and a = xy, yz, zx. 

(ii) For the (i 2 -site, we have 

|l,l,2,0) Q/3 = | t) Q | T)/3 

|l,0,2,0) Q/3 = -J=(| T)aU)/3 + U>«| T)/J) 

11,-1,2,0)^ = 1 l) a \l)p 
for 5=1, and 

|o,o J 2,o) 0/J = -^(| - U)«| T)/j) 

|0,0,2,l) a = |d) a 
for 5 = 0, where (i means the double occupancy. 

(iii) For the <i 3 -site, we have 
|3/2, 3/2, 3, 0) Q/ 3 7 = | ITT) 



11/2,-1/2,3,0)^= I I - ; I 1 , 2! | | |;,) 



|l/2,-l/2,3,l) a/ ,= |d) a | i)p 

for S = 1/2, where we use the notation \cr\O20D = 
1 01 ) Q I (72 ) /3 1 03 ) 7 . These states are used in the following as 
the intermediate states of the second-order perturbation 
processes. 

We here introduce an approximation; because the hop- 
ping parameters take the values t a ^> % ~ t c , 3 ' we 
assume t& = i c = for simplicity as in ref. 3 - 1 This means 
that the terms like H^ cx t a t c and fl^? cx tbt c are all 
neglected, retaining only the terms like Hf^ cx t\ in the 
second-order processes. We note that the orbital fluctu- 
ations are completely suppressed in this approximation 
because only two orbitals connected with the diagonal 
hopping t a comes out and no off-diagonal hopping terms 
appear in the effectve Hamiltonian. In other words, we 
obtain the effective spin-orbit Hamiltonian consisting of 
orbital-diagonal spin-sub-blocks with vanishing orbital 
off-diagonal blocks. 




Fig. 4. Schematic representation of five types of the bonds with 
different exchange interactions, (a) FM-1: the ferromagnetic 
bond with the process dfd 2 -> d|d] -» dfd 2 , (b) FM-2: the 
ferromagnetic bond with the process d 2 dh -t d'dj -> rf? , (c) 



3 1 J 

FM-3: the ferromagnetic bond with the process a~d - 



3 • 3 

d 2 d-, (d) AF-1: the antiferromagnetic bond with the process 
d 2 d 2 — > d?d}. or d 1 d 3 — » d 2 d 2 , and (e) AF-2: the antiferromag- 



netic bond with the process d 2 dh 



■ dftf or 



dU 2 

» 3 



d 2 d 1 . 

* 3 



|3/2, 1/2, 3, 0) Q/37 = -^(| TTI) + I TIT) + I iH> 



1 



13/2,-1/2, 3,0)^=^(1 UT) + UU) + I TU) 



13/2,-3/2,3,0)^ = I 444) 
for S = 3/2, and 

1 

72 
1 

75 

|1/2,1/2,3,1)^= \d) a \ Dp 



|l/2,l/2,3,0)g 7 =-5=(| nr )~ M, >> 

|i/2, 1/2,3,0)$, = -W| TIT) + 1 4TT) - 2| TTD) 



|l/2,-l/2,3,0)g 7 =-^(|T4l)-|ITI)) 
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We then have five types of bonds of spin exchange 
interactions as shown in Fig. 4; three of them are the 
bonds with ferromagnetic (FM) exchange interaction due 
to double-exchange mechanism and two of them are the 
bonds with antiferromagnetic (AF) exchange interaction 
due to kinetic-exchange mechanism. Defining the spin-1 
operator as S; and spin-1/2 operator as s,, we have the 
following Hamiltonian for each bond shown in Fig. 4. 

(i) The bond FM-1 with the process dfd 2 -> d^d) — > 



d\d): 



Elf = - JSi • S 3 + cl 



J = 



2f 



3(U'-V-J H ) 5(C/' -V + 2J H ) 
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2f 



2f 



5(U' -V + 2J H ) 



3(17' - V - J H ) 
where 1 is the unit operator. 



(ii) The bond FM-2 with the process djd) 



Wf = -2JS, 



cl 



J 



2f 



3(2C/' 



W - 
2^ 



2J ff ) 5(2C/' 



3V- 

2t? 



3(2Z7' - 3V - 2 J H ) 5(2Z7' - W + J# ) 
(iii) The bond FM-3 with the process dfdj 



did]: 



d\d) 



Hff = -2 JS, 



t 2 



cl 



3*| 
AV 



4(V + 2J H ) 



A(V + 2J H ) 

In the above three, we have the competing exchange 
couplings with positive and negative sign, but the sum 
is always positive for realistic values of the parameters 
J H , V, and U'. 



(iv) The bond AF-1 with the process dfd 2 



d]d* 



d\d): 



d\d] or 



H 



eff 



J Si " S A 



cl 



J =-c = 



U-V + Jh 
(v) The bond AF-2 with the process d\d) 



d\d) 



h: 



eff 



did}: 
= 2 JS 4 • s 7 - + cl 



d^d° or 



J = -c = 



2([7 + [/' - 37) 2(Z7 - [/' + V + Jh) 

These two exchange couplings have the doubly occupied 
states in their intermediate state and thus they are an- 
tiferromagnetic. 

Note that the intersite Coulomb repulsion V can be 
included in the bond Hamiltonian as shown above be- 
cause, due to the fixed CO pattern of Fig. 3, all the 
intermediate states have the same intersite Coulomb en- 
ergy, irrespective of the location of the bond. Also noted 
is that the so-called three-site terms in the perturbation 
do not appear for the CO pattern to be fixed. 

Thus, we have the effective spin-orbit Hamiltonian 



-ffeff = 



E 

<ij> 



off 



where the sum runs over all the nearest-neighbor pairs of 
sites. Note that this effective spin-orbit Hamiltonian has 
the form of block-diagonal in the spin<g>orbit space; i.e., 
orbital off-diagonal blocks are all zero. In other words, 
we have 3 N 00 patterns for the iV-site systems, and for 



each of them, we have the spin Hamiltonian. If we di- 
agonalize all the spin Hamiltonians, we can obtain the 
eigenstates of our effective spin-orbit Hamiltonian. We 
use the Lanczos diagonalization technique on small clus- 
ters to diagonalize the spin sub-block Hamiltonians and 
by comparing the lowest enrgies obtained we find the 
ground-state 00 patterns in the parameter space. 

Table I. Number of bonds in each OO pattern calculated for the 
12-site cluster. 



type of bond 


phase-I 


phase-II 


phase-Ill 


phase-IV 


FM-1 


4 


2 





6 


FM-2 


4 


2 





4 


FM-3 


12 


12 


12 


10 


AF-1 


4 


4 


4 


2 


AF-2 





2 


4 


2 




Fig. 5. (a) Calculated ground-state phase diagram of the effective 
spin-orbit Hamiltonian. The phases I— IV have the OO spatial 
patterns shown in (b). The region indicated as Stot > has the 
OO pattern I but has higher total spins. The CO of V 3 + and 
V 4+ is unstable in the upper-right part of the phase diagram, 
(b) Calculated OO spatial patterns of the phases I— IV. Shaded 
lobes have the electrons. 
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§5. Phase diagram 

Now, let us calculate the ground-state phase diagram 
of the effective Hamiltonian derived in the previous sec- 
tion. We use the 12-site (corresponding to 36-orbital) 
cluster with eight 5 = 1 spins and four S — 1/2 spins 
(corresponding to the filling of 20 electrons) coupled with 
the derived exchange interactions in the lattice of the CO 
pattern; we thereby calculated the ground state of each 
spin Hamiltonian. The periodic boundary condition is 
used. The calculated results for all possible 00 pat- 
terns (where we assume the unit cell containing 6 sites) 
indicate that a variety of the 00 patterns appear as 
the ground state, depending on the parameter values of 
Jh/U' and V/U' as shown in Fig. 5(a). 

The 00 spatial patterns of the phases I-IV are illus- 
trated in Fig. 5(b). Careful inspection of the patterns 
indicate that the gain in kinetic energy by the process 
FM-3 is the most important to stabilize these phases. 
This is evident in Table I, where the numbers of types of 
bonds existing in the phases I-IV are listed; i.e., the 00 
patterns I-IV are stabilized by maximizing the number 
of the bond FM-3 where the constant c-value is predom- 
inantly lower than others as shown in Fig. 6(a). 



1 ^R=8=g=2=2=8= s= = 8==2 H ! 

T * * • • " ■ i | 



i -2 



2 




FM-1 






-•- FM-2 






FM-3 


3 




-°~ AF-1 






~ A - AF-2 




0.40 



Fig. 6. Calculated vak,y^ pij^a^ ^jpyeg^change coupling constants 
J and (b) constant terms c in the effective spin Hamiltonians. 
The values normalized by the exchange coupling constant of the 
bond AF-1 are plotted. 



We find in Fig. 6(b) that the exchange energies of the 
antiferromagnetic bonds are generally larger than the ex- 
change energies of the ferromagnetic bonds; Jaf > Jf- 
This energy difference determines the detailed energy dif- 
ferences of the phases I-IV. For example, it may be that 
the strong AF-1 coupling promotes the singlet formation 
of the S — 1 spin pair by lowering its energy, to result 
in the appearance of the relatively stable phase I. Also, 
when the value of Jh is small, the strong antiferromag- 
netic couplings AF-1 and AF-2 stabilize the phase III 
as is noted in Table I. On the other hand, when Jh is 



large, the ferromagnetic couplings FM-1 and FM-2 sta- 
bilize the phase I as is also noted in Table I. When V/U' 
is small, the phase 5tot > 0, which has the same 00 
pattern as the phase I, becomes the ground state. This 
may be due to the strong ferromagnetic coupling FM-3, 
which is rapidly enhanced when V/U' becomes small as 
we find in Fig. 6(b). 

Because Jh/U' ~ 0.2 in the real material but the value 
of V/U' is not well-known, we may expect that the phases 
I-IV should all be possible to be realized. However, the 
experimental data on the spin degrees of freedom may 
single out a possible phase as we will discuss in the next 
section. 

§6. Discussion 

Finally, let us discuss the spin degrees of freedom of 
the system. Speculated spin states of the phases I-IV are 
illustrated in Fig. 7, which is supported in part by our 
small-cluster calculations of the spin correlation func- 
tions for the spin Hamiltonians. We find that the variety 
of the 00 patterns we have obtained result in a variety 
of the spin states as discussed below, although detailed 
numerical analyses of the Hamiltonians are needed for 
the definite descriptions of the states, which we leave for 
future study. 




Fig. 7. Speculated spin states of the phases I— IV, where shaded 
lobes have the electrons, dotted circles represent the local spin- 
singlet state, and arrows indicate the spin direction, (a) Phase- 
I: Singlets of two 5 = 1 spins and nearly-free 5 = 1/2 spins 
coexist, (b) Phase-II: Intermediate state between phases I and 
III. (c) Phase-Ill: State without spin frustration and thus with 
strong AF spin correlation, (d) Phase-IV: AF state with local 
high-spin clusters. 



For phase I, we find that the 00 spatial pattern may 
be consistent with the spin state where the spin structure 
may be of the partial singlet formation; i.e., 2/3 of V- 
ions have 5 — 1 and form the spin singlet pairs, leaving 
5=1/2 free spins on 1/3 of V-ions. For phase III, we 
find that the state has no spin frustration and thus we 
have the strong local antiferromagnetic spin correlations. 



6 



Yoshinori Shibata and Yukinori Ohta 



For phase II, we have the intermediate spin state between 
the phases I and III. For phase IV, we find that the state 
is formed by the local high-spin clusters of six V-ions, 
which arc linked in the antiferromagnetic arrangement. 

We may then suggest the following picture for the 
charge and orbital structure of the ground state of 
Bi x VgOi6: (i) There occurs the CO where the V-ions 
order as • • V 3+ V 3+ V 4+ - • • in the chain direction. Thus, 
one should observe the crystal structure of three-fold pe- 
riodicity, (ii) The possible spatial patterns of orbitals 
selected among the three i2g-orbitals are determined as 
the phase I in Fig. 5(b). (iii) The spin degrees of free- 
dom are then considered to be the partial singlet forma- 
tion. This state may be consistent with the result of an 
NMR experiment recently made by Waki et al.J' where 
it has been suggested that most of the spins form the 
singlet state with leaving a small amount of active spins, 
which then undergo an additional phase transition into 
the magnetic long-range order at lower temperatures. 

As for the MI phase transition of the real material, we 
suspect that the highly frustrated electronic state makes 
the system metallic at high temperatures, the CO occurs 
by lowering temperature to result in the MIT, and, trig- 
gered by this transition (but simultaneously) , there occur 
the 00 and spin-singlet formation. At lower temper- 
atures, there appears the magnetic long-range order, 2 ) 
which may be due to the weak exchange coupling be- 
tween the S — 1/2 spins in our 00 system (though not 
considered in this paper). 

The approximation used here neglects the orbital fluc- 
tuation completely The effect may however be impor- 
tant if we want to consider, e.g., the state at finite tem- 
peratures. The use of the approximation may however 
be justified because we here want to answer the question 
"what is the ground-state 00 pattern if we assume the 
presence of 00" . Possible absence of 00 as in LaTi03 8 ^ 
may of course be an interesting issue in the present sys- 
tem as well. Jahn- Teller distortions as a mechanism of 
the phase transition might also be relevant. Further ex- 
perimental (as well as theoretical) studies on the micro- 
scopic mechanism of this phase transition are therefore 
highly desirable. 

§7. Summary 

We have studied the electronic states of a vanadate 
material Bi^VsOig, a possibly new charge and orbital 
ordering system with the ^-orbitals on a ID triangular 
lattice with the mixed valency of 3d 2 : 3d 1 = 2 : 1 . By as- 
suming the charge ordering pattern, we have derived the 
effective spin-orbit Hamiltonian by the second-order per- 
turbation theory. Within the approximation neglecting 
small off-diagonal hopping parameters, we have found 
that the Hamiltonian is block-diagonal with vanishing 
orbital-off-diagonal sectors. We then have used a nu- 
merical diagonalization technique on small clusters of 
this Hamiltonian and have found that a variety of the 
ground-state phases with different orbital ordering pat- 
terns appear in the parameter space. We have argued 
that the orbital ordering pattern possibly realized in ex- 
periment should be a state of the partial singlets where 
the nearest-neighbor spin S = 1 pairs form singlets with 



leaving the nearly free S=l/2 spins. 

Although experimental data on this material are quite 
limited at present, we hope that the present theoreti- 
cal study will encourage further experimental studies to 
clarify the nature of the charge, orbital, and spin degrees 
of freedom of this interesting material. 
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